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Abstract 

m 

^ \ We point out that, if R-parity is broken spontaneously, the neutralino can decay 

to the final state majoron plus neutrino, which from the experimental point of view is 

. indistinguishable from the standard missing momentum signal of supersymmetry. We 

\Q \ identify the regions of parameter space where this decay mode is dominant and show 

^ ■ that they are independent of R-parity conserving SUSY parameters. Thus, (a) only 

. very weak limits on R-parity violating couplings can be derived from the observation 

\ of missing momentum events and (b) at future collider experiments huge statistics 

Q_i' might be necessary to establish that R-parity indeed is broken. Parameter combi- 

^ > nations which give calculated relic neutralino density larger than the measured dark 

D \ matter density in case of conserved R-parity are valid points in this scenario and their 

' phenomenology at the LHC deserves to be studied. 

: 

^ \ Is R-parity broken? Undoubtedly, a definite, affirmative answer to this question would 

^ \ have profound implications for both, particle physics and cosmology. Low-energy precision 
data can be used to derive limits on R-parity violating ) couplings [1] , but experimental 
proof that R-parity indeed is broken can only come from the observation of a decaying 
lightest supersymmetric particle (LSP) at a future collider experiment. 

In this paper we present and discuss a scenario with spontaneously broken R-parity where 
the model can be easily confused with the usual Minimal Supersymmetric Standard Model 
(MSSM) with conserved R-parity despite the fact that the LSP decays within a typical 
detector of present and future high energy experiments. We concentrate here on neutral 
LSPs as this is the case which has been studied mostly by experimental groups. However, 
in principle the LSP could be charged and quasi-stable on collider time scales and, thus, be 
observable. 

Suppose supersymmetric particles exist with masses, such that they can be produced 
at the LHC. To obtain evidence that the LSP decays, at least the following two conditions 
need to be met: (i) The decay length should not be (much) larger than the size of a typical 
detector. And, (ii) a minimum percentage of LSP decays, the exact number depending on 
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the SUSY spectrum and the luminosity of the colhder, of course, should lead to visible final 
states. It is straightforward to estimate that neutralino (x^) LSP decay lengths, for typical 
SUSY masses in the vicinity of (9(100) GeV, drop below L ~ 1 m for couplings larger 
than - very roughly - O(10~^). Condition (ii), however, is more subtle. Consider the lepton 
number violating part ^ of the explicit ^ superpotential [3] 

W = eiLiHu + XijkLiLjEl + X^-^LiQjDl (1) 

together with the corresponding bilinear soft terms, 

^oft = ^iBiLiHy^. (2) 

For trilinear , final states of x° decays always involve either two quarks and/or charged 
leptons. In bilinear ^p , scalar neutrino vacuum expectations values (vevs) induce couplings 
such as C^o^o,^ and Cy^y^ which lead to the decay x° 3u. Numerically [4], however, the 
branching ratio of this mode is always Br(x° 3u) < 0.1, since (a) Br{Z^ — > inv) ~ 0.2 
and (b) the coupling C^^o^o^ necessarily is C^ozo^ ~ Cy^O]Y±p. Thus, in explicit condition 
(ii) is automatically fulfilled. 

Spontaneous violation of R-parity (s-^p ) [5, 6], on the other hand, implies the existence 
of a Goldstone boson, usually called the majoron (J). In s--^pthe neutralino can then decay 
according to x° ~* J^) i-C- completely invisible. Below we will show that this decay mode 
can in fact easily be the dominant one, with branching ratios very close to 100 %. From the 
point of view of accelerator phenomenology, models with spontaneous violation of R-parity 
thus can resemble the MSSM with conserved R-parity and huge statistics might be necessary 
before it can be established that R-parity indeed is broken. 

Before turning to the details of the calculation, two more comments are in order. Firstly, 
experiments on solar and atmospheric neutrinos have shown that neutrinos have non-zero 
masses and mix [7]. Many models have been proposed, which potentially can explain these 
data. However, the presence of lepton number violation in ^p inevitably implies Majorana 
neutrino masses. Thus, significant upper bounds on most I/ip parameters follow from oscilla- 
tion experiments [8] as I/ip naturally implies a hierarchical neutrino mass spectrum. On the 
other hand, even the simplest bilinear I/ip model can explain measured neutrino masses [9], 
as can s-^p [10], if some specific combinations of I/ip parameters obey certain lower bounds, 
see below. In our numerical calculations we will always use ^p parameters fitted to neutrino 
physics. A short comment on how our results are affected by other choices of parameters 
is given near the end of this paper. 

Secondly, recent WMAP data [11] have confirmed the existence of non-baryonic dark 
matter and measured its contribution to the energy budget of the universe with unprece- 
dented accuracy. Because the neutrahno as a dark matter candidate is lost, most advocates 
of SUSY declare ^p undesirable. However, this does not mean that in models with ^p , 
the dark matter problem could not be solved. Possible DM candidates in I/ip are (i) light 

^Spontaneous IjLp conserves baryon number. One can also use baryon parity, or any similar symmetry, to 
get rid of the dangerous baryon number violating terms, see [2]. 
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gravitinos [12, 13], (ii) the axion [14, 15] or (iii) its superpartner, the axino [16], to mention 
a few. 

Turning to the definition of our model, we note that early versions of spontaneous Jjip 
used the left sneutrinos to break lepton number [5] . These models are now ruled out by LEP 
and astrophysical data [15, 17]. The model we consider [6] therefore extends the MSSM by 
three additional singlet superfields, namely, P^, S and with lepton number assignments 
of L = — 1, 1, respectively. The superpotential can be written as [6] 

+ KLi^H^-hQHaHu^ + h^v'S + ^^\ (3) 

The first three terms are the usual MSSM Yukawa terms. The terms coupling the lepton 
doublets to v'^ fix lepton number. The couphng of the field $ with the Higgs doublets 
generates an effective //-term a la Next to Minimal Supersymmetric Standard Model [18]. 
The last two terms, involving only singlet fields, give mass tov'^, S and once $ develops 
a vev. 

For simplicity we consider only one generation of d'^ and S. Adding more generations 
of or 5" does not add any quahtatively new features to the model. Note also, that the 
superpotential, Eq. (3), does not contain any terms with dimension of mass, thus potentially 
offering a solution to the //-problem of supersymmetry. The soft supersymmetry breaking 
terms of this model can be found in [10]. 

The basic guiding principle in the construction of Eq. (3) is that lepton number is con- 
served at the level of the superpotential. However, at low energy, i.e. after electro-weak 
symmetry breaking, various fields acquire vevs. Besides the usual MSSM Higgs boson vevs 
vd and vu, these are ($) = vp/\/2, (z/*^) = vr/\/2, (S) = t^s/V^ and (z/j) = vlJV^- Note, 
that Vji ^ generates efi^ective bilinear terms ej = h1vii/y/2 and that Vr, Vs and violate 
lepton number and R-parity spontaneously. 

Once R-parity is broken, Majorana neutrino masses are generated. After some algebraic 
manipulations one can express the effective (3,3) neutrino mass matrix in seesaw approxi- 
mation as [10] 

i'^^Z)ij = ^^i^j + H^i-^j + ^jK) + CeiEj , (4) 

where Aj = eiVd + //vl. and // = ^o^- The coefficients a, h and c can be found in [10]. 
Equation (4) resembles closely the structure found for the neutrino mass matrix in the 
explicit bilinear model at the one-loop level [9] . However, with the field content of Eq. (3) 
the spontaneous model produces two non-zero neutrino masses already at tree-level. Note, 
that the coefficients h and c go to zero (and thus one of the two neutrino masses vanishes at 
tree-level) in the limit when either h and ho vanish or, equivalently, vp ^ oo or vr ^ oo. 

Neutrino physics puts a number of constraints on the parameters Aj and e^. For the 
current paper the exact details are unimportant, however, the most essential constraint for 
the following discussion is that both, Aj/mgusv ^"^^ small. If some of the singlet 

fields are fight, i.e. have masses in the range of C(0.1 — few) TeV, |ei//t| can be as small as 
|ej///| ~ 10~^-10~^ (implying vlJvd ~ 10~^-10~^). On the other extreme, independent of 
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the singlet spectrum, |ej///| can not be larger than, say, |ej///| ~ 10~^, due to contributions 
from sbottom and stau loops to the neutrino mass matrix [9]. 

The requirement that f^. <^ Vd-, together with the assumption of vlJV < 1, can be used 
to find a rather simple expression for the majoron, in the basis A'q = Im{H^, i>i, i>2, 1^3, S, u'^): 

T I ~'"d'"L ^uVj Vli Vl2 Vl3 ^ Vs_ _Vr 

where = t>| + v\. We stress that neutrino data ensures vlJV < 1 even for vr as small 
as (9(0.1) GeV. Thus, the majoron in our model is mainly a gauge singlet. Constraints from 
the measurements of the invisible width of the boson [17] as well as all astrophysical 
bounds [15] are therefore easily satisfied. Note, that the majoron has no component in the 
$ direction. 

The coupling between the lightest neutralino, the neutrinos and the majoron in general is 
a complicated function of singlet and doublet parameters. However, in the limit vr, vs ^ 00 
one can derive a simple approximation formula, given by 

Cxlu^j ^ |iVu + '-^{g'Nn - gNu) + ■■■ (6) 

Here, e = C/J • e and vl — • vl, where C/J is the matrix which diagonalizes the neutrino 
mass matrix at tree-level. The dots in Eq. (6) stand for higher order terms. Equation(6) 
serves to show that for constant e and vl, C^q^^j as vr goes to infinity. This is as 
expected, since for vr —>■ 00 the spontaneous model approaches the explicit bilinear model. 
Note, that only the presence of the field P'^ is essential for the couphng Eq. (6). If S is absent, 
replace V ^ vr. 

We now turn to the numerical results. All calculations presented below are done with the 
program package SPheno [19], extended to include R-parity violation. We assume that the 
singlet states of the model are too heavy to be produced at future collider experiments such 
that only MSSM particles are observed. Note, that the smallness of the parameters guar- 
antees that R-parity violating decay modes of all supersymmetric particles can be neglected 
except for the LSP and that also SUSY particle production cross sections are completely 
MSSM-like. We will study only points where the lightest neutralino Xi is the LSP. 

As can be seen from Eq. (6), adjusting and vl. from neutrino data, see Eq. (4), one 
expects that the most important parameters for the decay x° — > Ju are {h"] = \/Si(^t)^ 
and, equivalently, vr. Therefore, in Fig. 1(a) we plot the sum over the branching ratios of 
neutralino decays containing at least one visible particle in the final state as a function of 
\h'^\, for fixed values of all other parameters except vp = 10-40 TeV. For this plot all MSSM 
parameters were chosen in such a way that the sparticle spectrum as well as all production 
cross sections very closely resemble the values of the SPA point SPSAla' [20]. For this point 
one expects about 2.5 10^ neutralinos at the LHC for a total luminosity of 300 fb~^. As 
Fig. 1(a) demonstrates, the total number of events with at least some visible energy in the 
final state drops in this example below 100 (10) for j/i'"! < 2.5 • 10~^ (7.5 ■ 10"^). Note, 
that in any SUSY event there are two Xi ^-iid that with such small branching ratios at most 
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Figure 1: (a) Sum over all at least partially visible decay modes of the lightest neutralino 
versus \n 

I = VT^iiKy for mSUGRA parameter point SPSla' [20]. (b) Same points, decay 
length versus Vji. This is the main result of the current paper. 



one of them decays visibly. Therefore, establishing that the Xi decays at the LHC becomes 
challenging. 

For completeness, in Fig. 1(b) we plot cr as a function of for the same points as in 
Fig. 1(a). For SPSla', the neutralino decays always with L < 1 cm. Other points have 
different upper limits, but L < 1 m always. 

We have explicitly checked that choosing different MSSM parameters reproduces Fig. 1(a) 
for a number of standard points from refs. [21, 22], as well as for points with a calculated 
fl^o ~ 1 (in case Rp were conserved). The invisible decay width of the neutralino does 
essentially not depend on MSSM parameters and thus is independent of ^dm obtained in 
conserved R-parity. 

For this model to resemble the MSSM it is not only important that the SUSY decays are 
MSSM-like and that the LSP decays invisibly but also that the Higgs sector behaves as in 
the MSSM. This is potentially in conflict with an interesting feature of this model, namely 
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Figure 2: Rji, defined as T{Sho JJ)/T{Sho — > bb) as a function of vr for the mSUGRA 
point SPSla' [20] for two values of the couphng h, for a number of different values of vp, 
vp = 10 TeV (blue) to vp = 40 TeV (red) in steps of 1 TeV. 

the invisible decay of the lightest doublet Higgs boson into two Majorons [10]. In Fig. 2 we 
show the ratio Rjb = T{Sho — > JJ)/T{Sho bb) for the same points as shown in Fig. 1, 
i.e. h = 0.009, as well as for a different choice of this parameter, h = 0.5. Note that the latter 
choice leads to invisible decay branching ratios indistinguishable from Fig. 1. Obviously 
there is enough freedom in choosing the parameters such that the invisible decay of the Higgs 
varies over several orders of magnitude, from a completely negligible level up to branching 
ratios similar to the dominant bb mode. We mention that the observed 'dips' in Rjb are due 
to a cancellation between various contributions in the Chojj coupling. Figure 2 demonstrates 
that not only the properties of the SUSY particles can be adjusted to be MSSM-like (with 
conserved R-parity) but also the Higgs sector. 

A final comment: We have chosen s-^p parameters fitted to neutrino masses. Smaller 
parameters e and vl lead to too small neutrino masses, but are allowed phenomenologically. 
These would lead to larger decay lengths, but not to larger branching ratios into visible final 
states 

In summary, we have studied the decay properties of a neutralino LSP in a model with 
spontaneous violation of R-parity. The neutralino decays invisibly in this model in large 
parts of the parameter space, such that even after several years of data taking at LHC it 
might be experimentally challenging to establish that R-parity is broken. A first hint of a 
non-standard scenario might come from cosmology, as in most of MSSM parameter space the 
relic density of the neutralinos is too large to be compatible with astronomical observations. 
However, with sufficiently large statistics it should eventually be possible to discover the 
rare visible decays of the lightest neutralino, e.g. into l^qq' final states, and thus prove that 
R-parity is violated. Needless to say, that such a search requires a dedicated effort by the 
experimentalists. 

^Adding some trilinear coupling, essentially unconstrained by neutrino masses [8], like A121 would of 
course lead to a larger number of visible events. 
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